A New Connection between Central Engine Weak Physics and the Dynamics of 

Gamma-Ray Burst Fireballs 
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We demonstrate a qualitatively new aspect of the dynamics of Gamma-Ray Burst (GRB) 
fireballs: the development of a substantial dispersion in the proton component in fireballs in 
■which neutron decoupling occurs and is sufficiently pronounced. This effect depends sensitively 
on the neutron to proton ratio in the fireball, becoming more dramatic with increasing neutron 
excess. Simple physical arguments and transport calculations indicate that the dispersion in 
Lorentz factor of the protons can be of order the final mean Lorentz factor of the fireball. 
We show how plasma instabilities could play an important role in the evolution of the fireball 
and how they might ultimately govern the development of such a velocity dispersion in the 
proton component. The role of these instabilities in setting/diminishing a proton Lorentz factor 
dispersion represents a new and potentially important venue for the study of plasma instabilities. 
Significant dispersion in the proton velocities translates into fewer protons attaining the highest 
Lorentz factors. This is tantamount to a reduction in the total energy required to attain a given 
Lorentz factor for the highest energy protons. As well, a proton component dispersion can have 
consequences for the electromagnetic and neutrino signature of GRBs. 

PACS numbers: 98.70.Rz, 14.60.Pq 



I. INTRODUCTION 

We sho-w that a simple, yet previously overlooked, 
mechanism may give rise to a substantial velocity dis- 
persion in the proton component of Gamma-Ray Burst 
fireballs. Namely, in GRB fireballs in -which neutron de- 
coupling occurs, large center of mass energy collisions 
bet-ween decoupled neutrons and protons will induce a 
velocity dispersion in the proton component. For certain 
fireball parameters the timescale for re-thermalization of 
these protons is longer than the dynamical timescale, and 
so the protons -will retain their acquired dispersion. 

This has consequences for the dynamics of the fireball. 
Since the ultra-relativistic proton shell is directly tied 
to gamma-ray photon production M 0| , the distribution 
of proton energy can affect the dynamics of photon emis- 
sion. In the internal shock GRB model the 7-ray photons 
are produced through internal multiple-shock collisions 
involving the protons. The efhciency for the conversion 
of shock kinetic energy to photon energy increases with 
increasing relative shock velocities [Q, g| . In addition, the 
energetic proton shock produces particle cascades lead- 
ing to multi-TeV neutrino emission y, g, Ig]. Therefore, 
both photon and neutrino radiation may be affected by 
the process discussed here. 

More basically, the process we discuss here represents 



a new aspect of our understanding of the evolution of 
homogeneous relativistic fireballs. The standard picture 
holds that a thermal plasma undergoes an initial ac- 
celeration phase, lasting until the energy in relativistic 
particles is transferred to the kinetic energy of protons, 
or, in the case of extremely baryon-dilute plasmas, until 
the plasma becomes optically thin to photons 0. Fol- 
lowing the acceleration phase a coasting phase ensues. 
In the coasting phase the baryons are assumed thermal 
with typically low temperatures, T < 10^^ MeV, in the 
plasma rest frame. 

Recently it has been recognized that in fireballs 
satisfying certain initial conditions the neutrons will 
dynamically decouple [g| |9|, |lO| from the expanding 
e /photon/proton plasma. There is an interesting con- 
nection between this decoupling and the net number of 
electrons per baryon. 
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where Y^. is the electron fraction, ng+ and n^- are the 
number densities of positrons and electrons respectively, 
and n and p are the number densities of neutrons and 
protons. The connection is that when conditions near the 
fireball source lead to a low Yg in the outflow, substantial 
differences between the final neutron and proton Lorentz 
factors are possible jlfl. We will demonstrate a second 
connection: when the electron fraction in the fireball is 
low, high energy collisions between decoupled neutrons 
and protons induce a velocity dispersion in the protons 
and that these "hot protons" are not re-thermalized. The 
sensitivity of this effect to Yg in the fireball is interesting 



because, in analogy to type II supernovae, the electron 
fraction may mirror weak physics in the central engine 

00. 



II. NEUTRON DECOUPLING IN 
RELATIVISTIC FIREBALLS 



The essential features of the physics of baryon flow 
in relativistic fireballs are obtained by considering a two 
component [(i) e /photon/proton and (ii) neutron] ho- 
mogeneous fireball with initial radius, temperature, and 
electron fraction Ro,To, and Yeo respectively. In terms of 
these quantities the entropy per baryon in units of lO^/cb 
is 



S5 « 1.250 X 10^2 ry(l MeV/To) 



(2) 



Here rj — Etot/M is the ratio of total energy to rest mass 
in the fireball. A central issue for GRBs is the "baryon 
loading problem" , the statement that rj must be large 
in order to get the fireball moving with a high Lorentz 
factor. 

For relativistic fireballs (?] > a few) , numerical and an- 
alytic work show that to first approximation the fireball 
evolves as ||l|, H 

7 = (To/T) = R/Ro for 7 < r; (3) 

7 = ?7 for i? > rjRo. (4) 

Here 7 is the Lorentz factor of the fireball. These re- 
lations follow from entropy and energy conservation and 
are violated at the beginning and end of the fireball's evo- 
lution (see below). In terms of the time t as measured 
by an observer comoving with the plasma, the Lorentz 
factor and temperature evolve as 7 = (Tq/T) — 6*/'^''='". 
Here Tdyn = Ro characterizes the fireball source size. (We 
adopt units where the speed of light is unity.) Obser- 
vations of time variability in GRBs give the constraint 
Tdyn < 10~^ s, while the smallest proposed GRB sources 
are solar mass scale compact objects with Tdyn > 10~^ s. 
Particles in the plasma suffer a 4-acceleration with 
magnitude gamma d'y/dR ~ Rq . Neutrons have such a 
small magnetic dipole moment that they are essentially 
only coupled to the accelerating plasma via strong scat- 
terings with protons (near the decoupling point anp ~ 
O.Ictt ^ lO^cTne ^ W^^an-y) whcrc ctt is the Thomson 
cross section. To an excellent approximation, then, the 
neutron-proton collision timcscalc in terms of the comov- 
ing proton number density np in the plasma rest frame 
is 



^p'^np'^rch 



(5) 



where Wroi is the relative neutron proton velocity. Apart 
from small corrections stemming from inelastic nucleon- 
nucleon scatterings, entropy is conserved throughout the 
acceleration of the fireball. This implies 



except for the brief period when the electron/positron 
pairs annihilate and transfer their entropy to the photons. 
Decoupling occurs if the dynamic and collision 
timescales become comparable before the end of the ac- 
celeration phase of the fireball's evolution. This condition 
is 
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(7) 



Here t^ is the dynamic timescale in units of 10~^ s, and 
(i;i.ciCTio)doc is the product of the relative velocity and 
neutron/proton cross section in units of 10 fm at the 
decoupling point. The precise v^-^i at which this latter 
quantity should be evaluated is unclear from our simple 
argument. This is not crucial because [vreicio] only varies 
by a factor of 5 from unity as Wrei increases from 10"'^ to 
1. The final Lorentz factor of the neutrons is estimated 
by evaluating 7 at the decoupling point. 
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Following decoupling, the proton/e /photon plasma 
begins to accelerate away from the neutrons. Energy 
conservation gives an estimate of the final mean Lorentz 
factor of the protons, 



(7p,final) ~ i^e ^ V' 7n,final 1 - Ye + 



^ E^ 7^ 
rrin 7n,final J 

(9) 

The term e(i?^/m„)(77r/7n, final) is present to account 
for energy lost to neutrinos arising from the decay of pi- 
ous created in inelastic nucleon-nucleon coUissions. Here 
TTin is the nucleon mass, e is the number of pious created 
per baryon (at most of order unity), Ei, is the average 
center of mass energy lost to neutrinos per pion decay, 
and 77r is the average Lorentz factor of created pious. 
For pious created in n-p collisions, the branching ratio 
for charged pion production below the two pion thresh- 
old is roughly 1/2, and E^ -- 50MeV. It is difficuh to 
analytically estimate "f,^, but as most inelastic nucleon- 
nucleon collisions occur near the decoupling point, it is 
presumably of order 7n, final- In the numerical examples 
we present below, neutrino energy loss is estimated and is 
found to be small enough so as not to alter the qualitative 
picture presented here. 



III. NEUTRON INDUCED PROTON HEATING 
IN RELATIVISTIC FIREBALLS 

Subsequent to neutron decoupling, protons will un- 
dergo collisions with decoupled neutrons. The first re- 
quirement for the presence of a "hot" proton component 
is that some of these collisions are high energy collisions. 
If in the lab frame a proton and neutron have Lorentz 
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factors 7p and 7n respectively, then the 7 of the neutron 
as seen by the proton is 



7rel = Jj(l + (7„/7pf) 

^ Tn 



(10) 



A high energy neutron-proton coUision iniphes then: 
(7p. final) /7n, final > afcw. This Condition is intcndend only 
as a rough guide for the fireball parameters needed to 
drive the proton component hot. 

For reasonable fireball parameters it is difficult to at- 
tain substantial (7p, final )/7n, final unless Ye < 0.5. It has 
been argued that a low electron fraction will naturally 
be obtained in many of the proposed GRB environments 
[ fiO| . In addition, in calculations of GRBs arising from 
neutron star mergers Yg « Q.l has been estimated [ [l4[ . 
It is not clear if such low electron fractions are also ob- 
tained, for example, in the coUapsar model for GRB cen- 
tral engines [ p^ . 

Protons undergoing collisions with neutrons will 
lose energy to the background photons and electrons. 
Whether or not these protons re-thermalize with the 
plasma depends on their thermalization timescale rthorm- 
The protons will remain "cool" {i.e., well coupled and 
thermal) if rtherm < Tdyn(< TcoIi), where the last inequal- 
ity holds near and after decoupling. If the opposite case 
holds, Tthcrm > ''dyn, the protons could keep the velocity 
dispersion they acquire through collisions with neutrons 
(i.e. become "hot"). 



mass (minus neutrino losses) goes into e"'"/e~ pairs and 
when, in addition, the processes creating e'^ /e~ pairs are 
sufhciently rapid that they come into equilibrium with 
electron/positron pair annihilation. 

By considering this limiting case, one finds that a 
pion-induced electromagnetic cascade can delay the on- 
set of a dispersion in the proton component (i.e. make 
Ttherm,pc {T = T'docoupio/4) < Tdyn) Only for vcry smaU 
electron fractions. Ye < 0.01. (The number 4 appear- 
ing in parentheses here arises from requiring that a large 
center of mass energy collision is needed to drive a proton 
dispersion and is explained in more detail below) . There- 
fore, a pion induced electromagnetic cascade has little or 
no leverage in delaying the development of a dispersion 
in the protons for moderate electron fractions. 

The second energy loss mechanism for protons arises 
because, when the protons acquire a dispersion, the elec- 
trons must rearrange themselves in order to preserve 
charge neutrality. This rearrangement, on its own, has 
a negligible effect on the proton dispersion because of 
the smallness of the ratio of the electron to proton mass. 
However, Compton scattering following this rearrange- 
ment can be a source of energy loss for the protons. In 
other words, runaway protons may be viewed as being 
tightly paired with electrons via the coulomb force, and 
the proton loses energy not only at the rate given in Eq. 
O, but in addition at the rate at which the electron to 
which it is coupled loses energy to the background pho- 
tons via Thomson/Compton drag: 



A. Collisional re-thermalization mechanisms 

There are two energy loss mechanisms for protons with 
a velocity relative to the background plasma. The first 
is proton-electron scattering. Bethe's formula |l6| gives 
this rate as: 
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Here ng is the electron number density in units of cm^"^ 
and V and 7 refer to the velocity and Lorentz factors 
respectively, of the proton relative to the plasma. 

The relations derived above imply that de- 
coupling occurs at a temperature Tdccoupic ~ 

TD(s5/(VrclfTlo)dccYeT5)l/3, with Tq = 0.005 McV. 

This clearly implies that decoupling occurs after e 
annihilation. Noting that the ionization electron num- 
ber density is Ue ~ 10^'^T^^yS'^^YeCm~^ , and noting 
that typical kinetic energies of protons relative to the 
plasma are of the order of a few GeV, we see that the 
thermalization timescale via proton/electron scattering 

is nhoonj,c « 8 ■ 10-^ S (Td/T)' SrjYe. 

Ref. made the observation that a pion-induced elec- 
tromagnetic cascade may lead to an increase in the num- 
ber density of electrons and positrons in the fireball. The 
greatest this increase could be occurs when all of the pion 
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with U-y the photon energy density, giving, rth« 
10^^s{Tj:)/T) . The total thermahzation timescale is 
given by 
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(13) 



A hot proton component arises, then, if neutron decou- 
pling occurs, and at the end of the fireball's acceleration 
phase Tthcrm > Tdyn- Applying the simple picture given 
above for the evolution of the fireball gives the second 
condition for a hot proton component: 



4>(4oo + y;)t-5. 



(14) 



In fact, the fireball is still accelerating for fireball radius 
values greater than rjRo (|12], and see below). This im- 
plies that the temperature during the acceleration phase 
becomes lower than Tq/tj and, hence, that the condition 
for the development of a hot proton component is some- 
what weaker than that given above. 

In Fig. 1 we illustrate the decrease in degree of ther- 
malization with increasing initial Tthcrm/''dyn for a proton 
moving relative to a background plasma and losin g en - 
ergy via the the processes described in Eqs. ^ and |lj. 
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FIG. 1: Illustration of the evolution of Lorentz factor of a 
test proton in a background plasma. The solid line is for 
the Lorentz factor of the background plasma and the dashed 
lines are for the Lorentz factors of the test protons. Different 
dashed lines correspond to different values of rthcrm/Tdyn (i.e. 
to different plasma temperatures) at lO^cm. The different 
values of rthcrm/Tdyn are given next to the curves. 



B. Collisionless considerations 

In addition to the requirement of charge neutrality 
leading to proton energy loss via electron Compton scat- 
tering, there are other collisionless processes that may af- 
fect the dynamics of the protons. Indeed, the proton and 
electron plasma frequencies are comparable to or larger 
than the collision frequencies after decoupling. 
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with Wpi and Wpo the proton and electron plasma fre- 
quencies. The post decoupling proton dispersed fireball 
may provide fertile and novel ground for the study of 
plasma instabilities. In particular, the protons may be 
subject to an analog of the two-stream instability be- 
cause of coUisionally produced bumps in the proton dis- 
tribution function at the bulk plasma and near the mean 
neutron component velocities. 

The plasma may also be subject to velocity space 
anisotropy-driven electromagnetic instabilities as the 
transverse velocity freezes out during the expansion of 
the fireball. Relevant for the present work, it has been 
shown that magnetic trapping can saturate these insta- 
bilities at substantial ion anisotropics |lj] , and that they 
are moderately stabilized by relativistic effects M, nol. 
A detailed study of the saturation level of the proton 
anisotropy in this system would be interesting because 
lower saturation levels imply less dispersion reduction 
due to adiabatic cooling of the tranverse momenta. In 
the presence of a magnetic field, either generated by the 
central engine or created through instabilities, other in- 
stabilities, involving, for example, ion-cyclotron waves, 
may operate. 



FIG. 2: A schematic picture of the various influences on the 
proton Lorentz factor distribution function during the accel- 
eration of the fireball. The solid line represents a coUisionally 
produced proton distribution function, while the dashed line 
represents the tendency of instabilities to smooth the distri- 
bution as collisionless processes liberate free energy. In this 
figure, pp and np represent proton-proton and neutron-proton 
collisions, while pe and pe7 represent the processes in Eqs. hil 
and |l^, respectively. 



Instabilities may play an important role in the set- 
ting of the proton distribution function. They are not 
expected to alter our result of a dispersion in the pro- 
ton component, because for the closely related system of 
counter-streaming ion beams they are well studied and a 
substantial reduction in an initial dispersion is not seen 
|17| , EG, El|. These instabilities have also been stud- 
ied in the regime where collisions are important, and in 
this case counter-streaming ion beams take, as expected, 
a few collision timescales to slow ||2^. However, stud- 
ies of the coupling between plasma instabilities and the 
hydrodynamic evolution of the protons in the peculiar 
high photon entropy/small dynamic timescale environ- 
ment characteristic of GRB fireballs have apparently not 
been done, and may give interesting and unexpected re- 
sults. The fact that plasma instabilities play a role dur- 
ing the acceleration phase of some relativistic fireballs is 
quite interesting, because the previous picture was that 
coUisional processes are sufficient to describe the early 
evolution in the absence of central engine-generated mag- 
netic fields. 

A diagram illustrating the various influences on the 
proton distribution function, including the trend toward 
a smoother distribution function resulting from plasma 
instabilities, is shown in Fig. 2. 



Dispersion in the proton component: analytic 
and numerical estimates 



When neutron decoupling occurs, and is pronounced, 
a rough estimate gives the dispersion in proton Lorentz 
factor as (7p, final) — 7n. final- (A somewhat more precise 
measure of dispersion will be given below.) For example, 
consider a fireball described by Y^ — 0.2, T5 = 2.7 and 
suppose that Tq — 3 and S5 = 4 (implying rj ^ 1000). 
We find 7„,finai ~ 340 and 7p, final ^ 3.4 x 10^ implying a 



dispersion in the proton component of 10 



V- 



We have defined the dispersion in terms of the central 
engine rest frame. For some purposes it is more useful to 
express the dispersion in terms of the mean proton rest 
frame. A simple approximate measure of the dispersion 
in this frame is 

1 / (TWmal) \ '/' (^„ ,„^^/(^p ,„^,) + 1) , (16) 

^ V 7n, final / 

which is at most of order a few. 

A more precise determination of the dispersion in the 
protons requires a transport calculation describing the 
evolution of the fireball. A simple single angular zone, hy- 
drodynamically consistent, steady state relativistic wind 
transport calculation has been performed. Neutrino en- 
ergy loss is estimated by assuming that following an 
inelastic collision the available kinetic energy is shared 
equally between two baryons and a pion. The pion en- 
ergy is then assumed lost from the system. In the calcu- 
lation, protons suffering large energy (7i.ci > 2) collisions 
with neutrons are assumed to decouple from the plasma 
if Ttherm > oTdyn , whcrc a IS taken as a free parameter. 
In the limit where neutrons and protons are coupled to 
the photon/e^ plasma our calculation reduces to that 
done in the pioneering study by Paczyhski ||23| on ultra- 
relativistic winds from compact sources. 

In Fig. 3 we present the results of a calculation for the 
parameters listed above: t^ = 2.7, Tq = 3, S5 = 4 and 
for the cases where Ye = 0.1,0.2,0.3, and 0.5. Results 
for different values of Yg have been presented to illustrate 
that the lower the electron fraction, the more pronounced 
the effects of collisions between neutrons and protons. 

To illustrate the effects discussed here we present re- 
sults for the case where (i) protons and neutrons are as- 
sumed coupled to and thermalized with the plasma; (ii) 
a distribution function is used to describe the neutrons 
but the protons are assumed coupled to and thermalized 
with the e /photon plasma (the agreement between the 
final proton Lorentz factors in this case and the simple 
estimate given by Eq. ^ is good); and (iii) a velocity 
distribution function is used to describe the neutron ve- 
locities, and protons are assumed to decouple from the 
plasma once they have undergone a collision with a neu- 
tron where 7rei > 2. For part (iii) of the calculation, 
the parameter a was taken to be 5. For the example 
parameters we have chosen the condition Ttherm > 5rdyn 
is always satisfied at the stage in the fireball's evolution 




FIG. 3: Evolution of proton Lorentz factor in a steady state 
relativistic wind characterized by the parameters given in the 
text, for different cases of electron fraction. Ye. The dot 
dashed line is for case (i), where baryons are assumed coupled 
to the plasma; the dashed line is for case (ii), where neutron 
velocities are described by a velocity distribution function and 
where neutrons are allowed to decouple while protons are as- 
sumed frozen into the plasma; and the solid line is for case 
(iii), where neutrons and protons are both allowed to decou- 
ple. Note the change in scale between the upper and lower 
figures. 



when 7i.ei > 2 collisions occur, and so the first condition 
is somewhat superfluous here. 

Perhaps the most striking feature of this calculation is 
the increase in Lorentz factor of the highest energy pro- 
tons as the protons acquire a dispersion. Because energy 
conservation implies that the mean Lorentz factor of the 
protons is unchanged between cases (ii) and (iii) (mod- 
ulo changes in the net energy carried by the neutrons 
and pions due to the protons acquiring a velocity disper- 
sion) the difference in final Lorentz factors between cases 
(ii) and (iii) provides a measure of the proton disper- 
sion. Note that because our calculation does not allow 
for a proton to interact with the expanding plasma once 
it has undergone a high energy collision with a neutron 
and once Ttherm > Srdyn our calculation may somewhat 
overestimate the proton dispersion. 

In Fig. 4 we display the evolution of a measure 
of dispersion in the proton component of the fireball 
(A7/7p_finai)- In this figure A7 is defined to be the differ- 
ence in proton Lorentz factors between models (ii) and 
(iii) for the case Ye = 0.2 shown in Fig. 1. 



IV. CONCLUSIONS 

We have pointed out and made estimates of an over- 
looked (and possibly generic) feature of the baryon (neu- 
tron/proton) flow in GRB fireballs: the acquisition of a 
large dispersion in Lorentz factor of the proton compo- 
nent subsequent to neutron decoupling. This represents 
a qualitatively new feature of GRB fireballs because the 
previous picture of the evolution of the proton component 




FIG. 4: Evolution of A7/{7p_finai) with radius for the case 
Ye = 0.2 and the parameters given in the text. 



held that the protons were thermal, even when neutron 
decoupling occurs. A dispersion in the proton component 
also implies that plasma instabilities may play an impor- 
tant role and suggests new studies of the role of collision- 
less processes during the acceleration and later shocking 
phase of the fireball's evolution. Because the effect we 
have described is principally sensitive to the neutron to 
proton ratio in the fireball, this work represents a new 
connection between the central engine weak physics that 
sets Ye and GRB fireball dynamics. Because shocks in- 
volving protons are thought to give rise to the photon 
signal characteristic of GRBs, a proton dispersion also 
has implications for the details and efficiency of photon 
production. A direct consequence for the charged parti- 
cle dynamics of neutron decoupling is interesting because 
the signal from the decay of pions created in inelastic col- 
lisions during neutron decoupling is expected to be weak. 
Detector event rates for the neutrinos from pion decay 
are estimated at a few per year p3. Direct evidence of 
neutron decoupling in a GRB event would provide in- 
formation about the progenitor fireball parameters and 
perhaps give a clue about the properties of the central 
engine. 
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